Relativistic Hartree-Bogoliubov model with density-dependent 

meson-nucleon couplings 

T. Niksic and D. Vretenar 
Physics Department, Faculty of Science, 
University of Zagreb, 10 000 Zagreb, Croatia 

g : P. Finelli 

O 



(N 



Physics Department, University of Bologna, 



>> and INFN - Bologna, I- AO 126 Bologna, Italy 

C3 



P. Ring 

Physik- Department der Technischen Universitat Miinchen, D-85748 Garching, Germany 

1 (Dated: February 8, 2008) 

O ' 

Abstract 

I The relativistic Hartree-Bogoliubov (RHB) model is extended to include density dependent 

meson-nucleon couplings. The effective Lagrangian is characterized by a phenomenological density 
dependence for the a, lo and p meson-nucleon vertex functions, adjusted to properties of nuclear 

73 

matter and finite nuclei. Pairing correlations are described by the pairing part of the finite range 



Gogny interaction. The new density-dependent effective interaction DD-ME1 is tested in the 



^ . analysis of the equations of state for symmetric and asymmetric nuclear matter, and of ground- 

h ; 

state properties of the Sn and Pb isotopic chains. Results of self-consistent RHB calculations are 
compared with experimental data, and with results previously obtained in the RHB model with 
non-linear self- interactions, as well as in the density dependent relativistic hadron field (DDRH) 
model. Parity-violating elastic electron scattering on Pb and Sn nuclei is calculated using a rel- 
ativistic optical model with inclusion of Coulomb distortion effects, and the resulting asymmetry 
parameters are related to the neutron ground-state density distributions. 

PACS numbers: 21.60.-n, 21.30.Fe, 21.65. +f, 21.10.-k 
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I. INTRODUCTION 



Models based on concepts of nonrenormalizable effective relativistic field theories and 
density functional theory provide a rich theoretical framework for studies of nuclear struc- 
ture phenomena, not only in nuclei along the valley of /3-stability, but also in exotic nuclei 
with extreme isospin values and close to the particle drip lines. A well known example of 
an effective theory of nuclear structure is Quantum Hadrodynamics (QHD), a field theo- 
retical framework of Lorentz-covariant, meson- nucleon or point-coupling models of nuclear 
dynamics M. The effective Lagrangians of QHD consist of known long-range interactions 
constrained by symmetries and a complete set of generic short-range interactions. The 
QHD framework implicitly includes vacuum effects, chiral symmetry, nucleon substructure, 
exchange terms, long- and short-range correlation effects. 

Structure models based on the relativistic mean-field (RMF) approximation have been 
successfully employed in studies of spherical and deformed nuclei all over the periodic ta- 
ble [|J. The relativistic framework has also been applied in studies of nuclear structure phe- 
nomena in nuclei far from (3- stability. In particular, in studies of isotopic chains that also 
included exotic nuclei with extreme isospin values, we have used the relativistic Hartree- 
Bogoliubov (RHB) model, which encloses a unified description of mean-field and pairing 
correlations. A number of interesting nuclear properties have been studied with the RHB 
model: the halo phenomenon in light nuclei || , properties of light neutron- rich nuclei M , 
the reduction of the effective single-nucleon spin-orbit potential in nuclei close to the drip- 
lines || , properties of neutron-rich Ni and Sn isotopes |J , the location of the proton drip- line 
between from Z = 31 to Z = 73 and the phenomenon of ground-state proton radioactivity 

The details of calculated nuclear properties that can be compared with empirical data, 
as well as the predictions of new phenomena far from /5-stability (halo nuclei, neutron 
skins, suppression of shell effects in neutron rich nuclei, ground-state proton radioactivity 
beyond the drip-line, the onset of exotic collective modes), crucially depend on the choice 
of the effective RMF Lagrangian in the ph channel, as well as on the treatment of pairing 
correlations. Several phenomenological parameterizations of the effective Lagrangian have 
been derived that provide a satisfactory description of nuclear properties along the /5-stability 
line, as well as in regions far from stability. These effective interactions are characterized by a 
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minimal set of model parameters: meson masses and meson-nucleon coupling constants. The 
most successful RMF effective interactions are purely phenomenological, with parameters 
adjusted to reproduce the nuclear matter equation of state and a set of global properties 
of spherical closed-shell nuclei. In most applications of the RHB model, in particular, we 
have used the NL3 effective interaction |TIJ for the RMF effective Lagrangian. Properties 
calculated with NL3 indicate that this is probably the best effective interaction so far, both 
for nuclei at and away from the line of /3-stability. 

The limitations of standard RMF effective interactions are, however, well known, even for 
nuclei close to the stability line. They are more pronounced in the isovector channel, which is 
poorly constrained by the available experimental data on ground-state properties of nuclei. 



For example, in a recent analysis of neutron radii in the framework of mean-field models [11 



it has been shown that conventional RMF models systematically overestimate the values of 
r n — r p- It is well known that they also predict an equation of state of neutron matter 
that is very different from the standard microscopic many-body neutron matter equation of 
state of Friedman and Pandharipande ||12|| . The parameterization of the isovector channel of 
an effective RMF interaction is, of course, extremely important for possible extrapolations 
to neutron-rich or proton-rich nuclei. For example, various non-relativistic and relativistic 
mean-field models differ significantly in the prediction of the two-neutron separation energies 
of the Sn isotopes with A > 132 |T^] and, therefore, of the exact location of the drip-line. 
There are also differences in the predicted location of the drip-line on the proton-rich side || . 
The properties of effective interactions are very important for the description of phenomena 
in nuclear astrophysics. The choice of the effective RMF Lagrangian directly determines the 
calculated properties of neutron stars (radii, surface crust) |14[ [151 , as well as the properties 
of nuclei that take part in the r-process or rp-process of nucleosynthesis. 

In order to overcome the limitations of standard RMF models, several solutions have 
been put forward. An obvious choice is to extend the minimal six-parameter RMF model 
by including additional interaction terms in the isoscalar, as well as in the isovector chan- 
nel. The complex many-body dynamics is effectively included in additional non-linear self- 
interactions. This approach has been studied both in the framework of meson-exchange 



models [TJ|, [Tj| [T7|, [L8j and with relativistic point-coupling models |19|, |2(J. Even though 
some interesting results have been obtained, especially in applications to nuclear astro- 
physics, the situation is not satisfactory. The problem is that the empirical data set of bulk 
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and single-particle properties of finite nuclei can only constrain six or seven parameters in 
the general expansion of an effective Lagrangian |21|]. One can, of course, include interaction 
terms that describe specific phenomena, but their coupling parameters and even their forms 
cannot be accurately determined in this way. In some cases even the signs of interaction 
terms are not determined by the data. The general expansion of an effective Lagrangian in 
powers of the fields and their derivatives can be controlled by the "naive dimensional analy- 
sis" (NDA) pi H, H, H. NDA tests the coefficients of the expansion for "naturalness", i.e. 
this method controls the magnitude of the coupling constants. Although NDA can exclude 
some interaction terms because their couplings would be "unnatural", it cannot determine 
the parameters of a model on a level of accuracy that is required for a quantitative analysis 
of nuclear structure data. 

Instead of including additional non-linear self-interaction terms in effective RMF La- 
grangians, another possibility is to formulate an effective hadron field theory with medium 
dependent meson-nucleon vertices. Such an approach retains the basic structure of QHD, but 
could be more directly related to the underlying microscopic description of nuclear interac- 
tions. In the density dependent relativistic hadron field (DDRH) model of Refs. EJ|, 26 



the medium dependence of the vertices is expressed by a functional of the baryon field 
operators. The meson-baryon coupling constants in nuclear matter are adjusted to the 
Dirac-Brueckner (DB) self-energies. A Lorentz-invariant functional is defined to project the 
nuclear matter results onto the meson-nucleon vertices of the effective DDRH model for fi- 
nite nuclei. In the early version [pJJ] of this model, the density dependence was only included 
in the field equations after the variation. However, a consistent treatment of medium effects 
implies a variation of the vertex functionals with respect to the baryon field operators, and 
this results in additional rearrangement self-energies in the single-nucleon Dirac equation 
|26fl . In Ref. the model was applied in the calculation of ground-state properties of 



doubly-closed shell nuclei. Density dependent a and uj meson couplings were used that were 
derived from DB calculations using the Bonn A, B, and C nucleon-nucleon potentials. It 
was shown that the inclusion of rearrangement self-energies is essential for a quantitative 
description of bulk properties and single-particle spectra. The model has been recently ex- 
tended to hypernuclei P?| , neutron star matter [E_HJ, and asymmetric nuclear matter and 
exotic nuclei . The density dependent interactions have been derived from the Groningen 
and Bonn-A nucleon-nucleon potentials. For finite nuclei, in particular, the quality of the 
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calculated properties are comparable with non-linear RMF models. 

In Ref. JJOJ Typel and Wolter introduced a phenomenological density dependence for the 
a, to and p meson-nucleon couplings, adjusted to properties of nuclear matter and some 
finite nuclei. The parameters of their DDRH model were also compared with coupling 
constants derived from DB calculations of nucleon self-energies. The model was used to 
study the equation of state of symmetric and asymmetric nuclear matter, and the ground 
state properties of semi-closed shell nuclei. Even though pairing correlations were only 
treated in the BCS approximation, properties at the proton and neutron drip-lines were 
calculated. It was, however, emphasized that a more realistic description of nuclei at the 
drip-lines necessitates the use of the Hartree-Bogoliubov framework with finite range pairing 



interactions. The phenomenological ansatz of Ref. [|30| for the functional form of the density 
dependence of the meson-nucleon vertices was also used in Ref. |2{| to derive the meson- 
nucleon coupling parameters from the Groningen and Bonn-A NN potentials. 

In this work we present an extension of the relativistic Hartree-Bogoliubov (RHB) model 
that includes density dependent meson-nucleon couplings. For the effective RMF Lagrangian 
we follow the approach of Typel and Wolter [BO] and use their phenomenological functional 
forms for the density dependence of the vertex functions. The parameters of the effective 
interaction are, however, adjusted in a different way. Pairing correlations are described by the 
pairing part of the finite range Gogny interaction. The new model is tested in the analysis 
of the equation of state for symmetric and asymmetric nuclear matter, and ground-state 
properties of the Sn and Pb isotopic chains. The results are compared with experimental 
data and with previous results obtained in the RHB, as well as DDRH frameworks. 

In Sec. II we outline the RHB model with density dependent meson-nucleon couplings. 
The new parameterization of the vector density dependence of the vertex functionals is 
discussed in Sec. |H| in comparison with the effective interaction of Typel and Wolter, and 
we also analyze the results for the equation of state for symmetric and asymmetric nuclear 
matter. In Sec. |V| the new density-dependent effective interaction DD-ME1 is employed 
in the RHB calculations of ground-state properties of Sn and Pb nuclei. Binding energies, 
charge radii, differences between neutron and proton radii, spin-orbit splittings, and charge 
isotope shifts are compared with available experimental data, and with the results obtained 
with the non-linear interaction NL3. Sec. |V| contains an analysis of parity-violating elastic 
electron scattering on 208 Pb and on those Sn isotopes for which there are experimental data 
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on r n — r p . For the elastic scattering of 850 MeV electrons on these nuclei, the calculated 
parity-violating asymmetry parameters are related to the Fourier transforms of the neutron 
density distributions. In Sec. [VT| we summarize the results of this work and present an 
outlook for future applications of the RHB model with density-dependent meson-nucleon 
couplings. 

II. RELATIVISTIC HARTREE-BOGOLIUBOV MODEL WITH DENSITY- 
DEPENDENT MESON-NUCLEON COUPLINGS 

The framework of density- dependent hadron field theory is described in great detail in 
Refs. . In this section we outline the essential features of the model with vector 

density dependence of the meson-nucleon couplings. The model is defined by the relativistic 
Lagrangian density 

C = ip (ry • d — m) ip + -(da) 2 — -m^a 2 

- -n^QT + -m 2 u 2 - -R^BT + -m 2 p 2 - -F ul/ F^ 

- g a ^aip - ■ wip - g p ^ ■ prip - e^7 • A ^ 1 Tp- ^ ■ (!) 

Vectors in isospin space are denoted by arrows, and bold-faced symbols will indicate vectors 
in ordinary three-dimensional space. The Dirac spinor ip denotes the nucleon with mass m. 
m a , rriu, and m p are the masses of the cx-meson, the u;-meson, and the p-meson. g a , g u , and g p 
are the corresponding coupling constants for the mesons to the nucleon. e 2 /47r = 1/137.036. 
The coupling constants and unknown meson masses are parameters, adjusted to reproduce 
nuclear matter properties and ground-state properties of finite nuclei, fl^", R^ v , and F^ v 
are the field tensors of the vector fields u, p, and of the photon: 

QT = d^u" - cTa* (2) 
BT = d»p v - d u p» (3) 
= d"A v - d u A» . (4) 

go, g^ji and g p are assumed to be vertex functions of Lorentz-scalar bilinear forms of the 
nucleon operators. In most applications of the density-dependent hadron field theory the 
meson-nucleon couplings are functions of the vector density 

Pv = with j p = ^7 M V • (5) 
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Another obvious choice is the dependence on the scalar density p s = ipip. It has been shown, 
however, that the vector density dependence produces better results for finite nuclei ||26|| , and 



provides a more natural relation between the self-energies of the density- dependent hadron 
field theory and the Dirac-Brueckner microscopic self-energies |29|| . In the following we 
assume the vector density dependence of the meson-nucleon couplings. The single-nucleon 
Dirac equation is derived by variation of the Lagrangian ([I]) with respect to ip 

[y*(i^-S M )-(m-E)]V = 0, (6) 

with the nucleon self-energies defined by the following relations 

£ = 9*° (7) 

(1 — T ) 

S M = guUp + g P r ■ p M + e y ^ 3) A p + E J . (8) 

The density dependence of the vertex functions g a , g u , and g p produces the rearrangement 
contribution £^ to the vector self-energy 

K = 3 - (^Wv^ + %*w*i> ■ Pu + . (9) 

M Pv \op v dp v dp v J 

The inclusion of the rearrangement self-energies is essential for the energy-momentum con- 
servation and the thermodynamical consistency of the model [ PU| , 50 ]. 

The lowest order of the quantum field theory is the mean-field approximation: the meson 
field operators are replaced by their expectation values. The A nucleons, described by a 
Slater determinant |$) of single-particle spinors ipi, (i = 1,2, A), move independently 
in the classical meson fields. The sources of the meson fields are defined by the nucleon 
densities and currents. The ground state of a nucleus is described by the stationary self- 
consistent solution of the coupled system of Dirac and Klein-Gordon equations. Due to 
charge conservation, only the 3-component of the isovector rho meson contributes. For an 
even-even system the spatial vector components oj, ps and A vanish, and the self-energies 
are determined by the solutions of the Klein-Gordon and Poisson equations 

(-A + m«Xr) = -g aPs {r), (10) 

(-A + m w )w(r) = g^Pvir), (11) 

(-A + m p )p 3 (r) = g p (p n (r) - p p (r)), (12) 

-Mo(r) = e 2 p c (r) . (13) 



The source terms on the left-hand side of these equations are sums of bilinear products 
of baryon amplitudes. The densities are calculated in the no-sea approximation, i.e. only 
occupied single-nucleon states with positive energy explicitly contribute to the nucleon self- 
energies. 

In addition to the self-consistent mean-field potential, pairing correlations have to be 
included in order to describe ground-state properties of open-shell nuclei. For spherical and 
deformed nuclei not too far from the stability line, pairing is often treated phenomenologi- 
cally in the simple BCS approximation 0. However, the BCS model presents only a poor 
approximation for exotic nuclei far from the valley of /^-stability. The structure of weakly 
bound nuclei necessitates a unified and self-consistent treatment of mean-field and pairing 
correlations. In particular, the relativistic Hartree-Bogoliubov (RHB) model p, |], || H 
represents a relativistic extension of the Hartree-Fock-Bogoliubov (HFB) framework. In the 
RHB model the ground state of a nucleus |$ > is represented by the product of independent 
single-quasiparticle states. These states are eigenvectors of the generalized single-nucleon 
Hamiltonian that contains two average potentials: the self-consistent mean-field T, which en- 
closes all the long range particle-hole (ph) correlations, and a pairing field A, which sums up 
the particle-particle (pp) correlations. In the Hartree approximation for the self-consistent 
mean field, the relativistic Hartree-Bogoliubov equations read 

tf*r\ v \'l (14) 

V V h (r) ) 

where ho is the single-nucleon Dirac Hamiltonian, and m is the nucleon mass. The chemical 
potential A has to be determined by the particle number subsidiary condition, in order that 
the expectation value of the particle number operator in the ground state equals the number 
of nucleons. The column vectors denote the quasiparticle wave functions, and Ej, are the 
quasiparticle energies. The source terms in equations (|10|) to (|1^) are sums of bilinear 
products of baryon amplitudes 

Ps(r) = E ^W7V fc (r), (15) 

E k >0 

Pv{r) = E VfrWr), (16) 

E k >0 

Pn (r)- Pp (v) = £ ^(r)r 3 y,(r), (17) 




E k >0 
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Pc(r) = £ ^(r)i-^(r), (18) 

£ fe >o z 

where the sums run over all positive energy states. In most applications of the RHB model 
a phenomenological pairing interaction has been used, the pairing part of the Gogny force, 



V pp (l,2) = e- ((ri ~ r2)/ ^ )2 (Wi + BiP a -HiP T -M i P a P r ), (19) 

i=l,2 



with the set D1S |31|] for the parameters fa, Wi, Bi, Hi, and Mi (i = 1,2). 

The RHB equations are solved self-consistently, with potentials determined in the mean- 
field approximation from solutions of Klein-Gordon equations for the meson fields. The 
Dirac-Hartree-Bogoliubov equations and the equations for the meson fields are solved by 
expanding the nucleon spinors Z7jt(r) and V&(r), and the meson fields in terms of the eigen- 
functions of a spherical or deformed axially symmetric oscillator potential. A detailed de- 
scription of the relativistic Hartree-Bogoliubov model for spherical and deformed nuclei can 
be found in Refs. [|J and ||, respectively. 

III. PARAMETERIZATION OF THE DENSITY DEPENDENCE OF THE 
MESON-NUCLEON COUPLINGS 

The density dependence of the meson-nucleon couplings can be obtained from microscopic 
Dirac-Brueckner (DB) calculations of nucleon-self energies in symmetric and asymmetric 



nuclear matter |26], |32|. However, depending on the choice of the nucleon-nucleon potential 
and on the approximations in the DB calculation, rather different results are obtained for the 
density dependence of the vertex functions. At low densities, in particular, DB calculations 
of nuclear matter become unreliable, and meson-nucleon couplings determined directly from 
DB self-energies provide only a qualitative description of ground-state properties for finite 
nuclei. Instead of adjusting the vertex functions directly to DB self-energies, in Ref. ]3(J an 
ansatz was made for the functional form of the density dependence that encloses different DB 
results. The parameters of the density dependence were obtained from a fit to properties 
of nuclear matter and finite nuclei. The same functional form was also used in Ref. |29| 



to adjust the meson-nucleon couplings to the DB self-energies derived from the Groningen 
and Bonn-A nucleon-nucleon potentials. In this work we adopt for the density dependence 
of the meson-nucleon couplings the functionals of Ref. [jSPR , and adjust the parameters to 
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properties of symmetric and asymmetric nuclear matter, binding energies, charge radii and 
neutron radii of spherical nuclei. The coupling of the a-meson and a;-meson to the nucleon 
field reads 

9i(p) = 9i(p sa t)fi(x) for % = a,u , (20) 

where 



is a function of x = p/p sa t- The eight real parameters in (|21[) are not independent. The 
five constraints fi(l) = 1, f"(X) = and /f (0) = 0> reduce the number of independent 

parameters to three. Two additional parameters in the isoscalar channel are g a {p S at) and 
^(Psat)- The functional form of the density dependence of the p- meson coupling is suggested 
by DB calculations of asymmetric nuclear matter [02] 



9 p{p) = tfpOsat) exp [-a p (x - 1)] . (22) 
The isovector channel is parameterized by g p (p s &t) and a p . 



In Ref. [[!(J the standard free values for the masses of the u and p mesons were taken: 
m uj — 783 MeV and m p = 763 MeV. The mass of the a meson was fixed to m a = 550 
MeV, and the remaining seven independent parameters were adjusted to nuclear matter 
properties (three parameters) and to the binding energies of symmetric and neutron-rich 
nuclei (four parameters). In this work the density-dependent meson- nucleon couplings have 
been determined using a somewhat different procedure. The parameters have been adjusted 
simultaneously to properties of nuclear matter (see Table |), and to binding energies, charge 
radii and differences between neutron and proton radii of spherical nuclei (see Table pl|) . In 



addition to the seven coupling parameters, the mass of the cr-meson has also been included 
in the fitting procedure, i.e. we have used one more free parameter with respect to the 



model of Ref. [30]. For the open shell nuclei in Table ID, pairing correlations have been 
treated in the BCS approximation with empirical pairing gaps (five-point formula). For 
nuclear matter the "empirical" input was: E/A = —16 MeV (5%), po = 0.153 fm~ 3 (10%), 
Kq = 250 MeV (10%), and J = 33 MeV (10%). The values in parentheses correspond to 
the error bars used in the fitting procedure. The binding energies of finite nuclei and the 
charge radii were taken within an accuracy of 0.1% and 0.2%, respectively. Due to large 
experimental uncertainties, however, the error bar used for the neutron skin was 5%. After 



10 



the solution of the self-consistent equations, we subtract from the total binding energy the 
microscopic estimate for the center-of-mass correction 

< P L > 

Ecm = 2A^T ' (23) 

where P cm is the total momentum of a nucleus with A nucleons. The resulting parameters 
of the density-dependent meson-exchange effective interaction (DD-ME1) are displayed in 
Table |, in comparison to those of Ref. |30| (TW-99). We note two additional differences 
between these models. In the relativistic mean-field model of Ref. |30| the center-of-mass 



correction (23) is calculated by using the non-relativist ic approximation for the nucleon wave 
functions, and the Coulomb energy is corrected by multiplying the vector self-energy of the 
proton with the factor (Z — 1)/Z. In Table | we first note that, even though m a is a free 
parameter in our model, the adjustment to nuclear matter and to properties of finite spherical 
nuclei produces a value that is very close to that of the TW-99 parameterization. The two 
effective interactions display similar values for the meson-nucleon coupling parameters at 
saturation density (?i(p S at) (i — cr ) ^ ) p) ) a s well as the value of the parameter a p , which 
determines the density-dependence of the p- meson coupling fl2"2"|). The eight parameters, 
which characterize the density dependence ([H]) of the a- and u;-meson couplings in DD- 
ME1, are very different from those of the TW-99 effective interaction. The reason is, of 
course, that only three of these parameters are independent. In our model b a , d a and 

are adjusted to properties of nuclear matter and finite nuclei. Nevertheless, the overall 
density dependence of the meson-nucleon vertex functions is very similar in the two models, 
as shown in Fig. [3] for the relevant densities p < 0.3 fm -3 . 

Nuclear matter properties calculated with the DD-ME1 interaction are illustrated in 
Table [II] and Figs. |2|, |3], and |4j. The results are shown in comparison with those obtained 



with the density-dependent effective interaction TW-99 | 30fl , and with two standard non- 



linear parameter sets NL3 |10| and NL1 | 3~3fl . The later non-linear effective interactions have 



been used extensively in studies of nuclear structure phenomena over the whole periodic 
table, from light nuclei to superheavy elements. For symmetric nuclear matter all four 
interactions display similar saturation densities (with NL3 at the low end), and binding 
energies per nucleon (with NL1 at the high end). While three parameters of TW-99 have 
been specifically adjusted to the values of p sat , E/A and the incompressibility K shown 
in Table O, the effective interactions DD-ME1, NL3 and NL1 have been simultaneously 
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adjusted to properties of nuclear matter and finite nuclei. The incompressibility modulus 
Kq of the two density-dependent interactions (~ 240 — 245) MeV lies between the values 
predicted by the non- linear interactions NL1 and NL3. Calculations of the excitation energies 
of isoscalar giant monopole resonances in spherical nuclei in the time-dependent relativistic 
mean-field framework and in the relativistic random-phase approximation J35], suggest 



that the nuclear matter incompressibility modulus should be in the range Kq « 250 — 270 
MeV. All four effective interactions display low values of the effective mass m* that are, of 
course, necessary in order to reproduce the empirical spin-orbit splittings in spherical nuclei. 
The equations of state of symmetric nuclear matter are compared in Fig. |2|. All four binding 
energy per particle curves display a very similar dependence on density below the saturation 
point p sat . Pronounced differences show up at higher densities. In particular, the two non- 
linear effective interactions NL1 and NL3 display a much steeper increase of the binding 
energy, especially with respect to TW-99. This is due to the fact that the cu-meson, which 
dominates at higher densities, enters linearly in NL1 and NL3, with constant coupling. The 
equation of state calculated with the DD-ME1 interaction shows an intermediate density 
dependence in this region, though closer to NL1 and NL3. 

The principal difference between the density-dependent effective interactions DD-ME1 
and TW-99 on one hand, and the non-linear interactions NL3 and NL1 on the other, are the 
properties of asymmetric matter. This is a very important point, because different isovector 
properties in nuclear matter lead to very different predictions for the properties of exotic 
nuclei with extreme isospin values. The energy per particle of asymmetric nuclear matter 
can be expanded about the equilibrium density p sat in a Taylor series in p and a [[36 



where 



E(p, a) = E(p, 0) + S 2 (p)a 2 + S 4 (p)a 4 + ■■■ (24) 



N-Z 



£(p,0) = -o„ + T |^- (p - p„) 2 + ... (26) 

-'sa 

Psat 18 Psat 



and 



ISpL 

5*2 (p) = a 4 + (p - Psat) + ft K 2 ° (p - Psat) 2 H (27) 
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The empirical value at saturation density S^Asat) = «4 = 30 ± 4 MeV. The parameter po 
defines the linear density dependence of the symmetry energy, and Ai^o is the correction to 
the incompressibility. The contribution of the term Si{p)a A in (|2^) is very small in ordinary 
nuclei and the coefficient is not constrained in the mean-field approximation. We first note 
that the non-linear effective interactions NL1 and NL3 have a considerably larger value a± 
of the symmetry energy at saturation density. This is also true for other standard non-linear 
parameter sets, and is due to the fact that the isovector channel of these effective forces is 
parameterized by a single constant, the density-independent p-meson coupling g p . With a 
single parameter in the isovector channel it is not possible to reproduce simultaneously the 
empirical value of 04 and the masses of N 7^ Z nuclei. This only becomes possible if a density 
dependence is included in the p-meson coupling, as it is done in TW-99 and DD-ME1. In 



a recent analysis of neutron radii in non-relativistic and covariant mean- field models [ I I 



Furnstahl has studied the linear correlation between the neutron skin and the symmetry 
energy. In particular, he has shown that there is a very strong linear correlation between 
the neutron skin thickness in 208 Pb and the individual parameters, which determine the 
symmetry energy S^p): 0,4, po and AK$. The empirical value of r n — r p in 208 Pb (0.20 ±0.04 
fm from proton scattering data [[IT], and 0.19 ± 0.09 fm from the alpha scattering excitation 
of the isovector giant dipole resonance ||5B| ) places the following constraints on the values of 
the parameters of the symmetry energy: a 4 « 30 — 34 MeV, 2 Mev/fm 3 < p < 4 Mev/fm 3 , 
and —200 MeV < AKq < —50 MeV. In Table [Tl] we notice that, while these constraints 
are satisfied by the density-dependent interactions DD-ME1 and TW-99, the parameters 
of the symmetry energy of the non-linear interactions are systematically much larger. In 
particular, p is too large by a factor « 2, and the correction to the incompressibility 
AKq has even a wrong sign for the two non-linear interactions. The qualitatively different 
density dependence of the symmetry energy for the two classes of effective interactions 
is also illustrated in Fig. ^ where we plot the coefficient S2 as a function of the baryon 
density. Due to the very large value of po and the small absolute value of AK , for NL3 and 
NL1 S2 displays an almost linear density dependence of p. For the two density-dependent 
interactions, on the other hand, the quadratic term of 5*2 dominates, especially at densities 
p > 0.1 fm~ 3 . 

In Fig. ^ we display the energy per particle of neutron matter as a function of the 
neutron density. At low densities, which are relevant for nuclear structure problems, the 
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results for the four relativistic mean-field interactions DD-ME1, TW-99, NL3 and NL1 are 
shown in comparison with the microscopic many-body neutron matter equation of state of 
Friedman and Pandharipande []I2| . The later is well reproduced by the density-dependent 
effective interactions, especially by TW-99, while the two non-linear interactions NL3 and 
NL1 display a qualitatively different neutron matter equation of state, even at very low 
densities. 

The density-dependent effective interaction DD-ME1 has been simultaneously adjusted 
to nuclear matter properties and to ground-state properties of the spherical nuclei shown in 



Table |IH| . The calculated binding energies, charge radii and differences between neutron and 
proton radii are shown in comparison with available experimental data. The choice of the 
isotopes of Sn and Pb for the fitting procedure was motivated by the desire to construct an 
effective interaction that could be applied in the description of long isotopic chains, including 
exotic nuclei which lie very far from the valley of /3-stability. The overall agreement between 
the calculated quantities and experimental data in Table |T| is very good. 

One of the advantages of using the relativistic framework lies in the fact that the effec- 
tive single-nucleon spin-orbit potential arises naturally from the Dirac-Lorentz structure of 
the effective Lagrangian. The single-nucleon Hamiltonian does not contain any adjustable 
parameter for the spin-orbit interaction. In Table |I^ we compare the energy spacings of 
spin-orbit partners in the doubly closed-shell nuclei 16 0, 40 Ca, 48 Ca, 132 Sn and 208 Pb, with 
the values calculated with the DD-ME1 interaction and with the prediction of the standard 
NL3 non-linear interaction. The experimental data are from Ref. fl39|l . We notice that, even 
though the values calculated with NL3 are already in very good agreement with experimen- 
tal data, a further improvement is obtained with the DD-ME1 interaction, especially for the 
lighter nuclei 16 0, 40 Ca and 48 Ca. 

IV. GROUND STATES OF THE Sn AND Pb ISOTOPES 

In Ref. |5| we have applied the relativistic Hartree Bogoliubov (RHB) model in a detailed 



analysis of ground-state properties of Ni and Sn isotopes. The NL3 parameter set flO 



was used for the effective mean-field Lagrangian, and pairing correlations were described 



by the pairing part of the finite range Gogny interaction D1S pi ]. Fully self-consistent 



RHB solutions were calculated for the isotopic chains of Ni (28 < N < 50) and Sn (50 < 



14 



iV < 82). Binding energies, neutron separation energies, and proton and neutron rms radii 
were compared with experimental data. The reduction of the spin-orbit potential with the 
increase of the number of neutrons was studied, and the resulting energy spacings between 
spin-orbit partners were discussed, as well as pairing properties calculated with a finite range 
effective interaction in the pp channel. 

In this section we test the new density- dependent meson-exchange effective force DD-ME1 
in comparison with the non-linear interaction NL3. The RHB model is used to calculate 
ground-state properties of Sn and Pb isotopes. Both NL3 and DD-ME1 mean-field La- 
grangians are employed for the ph channel, and the pairing part of the Gogny interaction 
D1S is used in the pp channel. This pairing interaction is a sum of two Gaussians with 
finite range and properly chosen spin and isospin dependence. The Gogny force has been 
very carefully adjusted to the pairing properties of finite nuclei all over the periodic ta- 
ble. Its basic advantage is the finite range, which automatically guarantees a proper cut-off 
in momentum space. By comparing results of fully self-consistent RHB calculations with 
experimental data, we will show that the new effective interaction DD-ME1 provides an 
excellent description of ground-state properties and, as compared with NL3, the isovector 
channel is considerably improved. 

In Fig. [| we plot the deviations of the theoretical masses of Sn isotopes, calculated in the 
RHB model with the DD-ME1 and NL3 interactions, from the empirical values [I0|. Both 
interactions display very good results over the entire major shell 50 < iV < 82. For the new 
interaction DD-ME1, in particular, only in few cases the absolute deviation of the calculated 
mass exceeds 0.1%. 

The isotopic dependence of the difference between the theoretical and experimental charge 



radii |4T| of Sn nuclei is displayed in Fig. y. The charge radii calculated with both DD-ME1 
and NL3 interactions are systematically smaller than the experimental values. The new 
density-dependent force, however, reduces the deviations from the experimental radii by a 
factor « 2. The parameters of DD-ME1 have been adjusted to the charge radii of 112,116 ' 124 Sn, 



and the absolute deviations for these nuclei can be compared in Table pi . 

The calculated differences between radii of neutron and proton ground-state distributions 
of Sn nuclei are shown in Fig. [7]. The non-linear interaction NL3 systematically predicts 
larger values of r n — r p . This effect is even more pronounced for the older parameter set 



NL1 IE2I. The difference between the values calculated with NL3 and DD-ME1 increases 
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with the number of neutrons to about 0.1 f m at iV = 82, but then it remains practically 
constant for iV > 82. The calculated values of r n — r p are compared with experimental 
data H3j in Fig. H. While both interactions reproduce the isotopic trend of the experimental 



data, NL3 obviously overestimates the neutron skin. The values calculated with DD-ME1, 
on the other hand, are in excellent agreement with the experimental data. This result 
presents a strong indication that the isovector channel of the effective interaction DD-ME1 
is correctly parameterized. 

In Refs. H, P it has been shown that the relativistic mean-field framework predicts a 
strong reduction of the magnitude of the spin-orbit term in the effective single nucleon 
potential of nuclei with extreme isospin values. Starting from T z = nuclei, and increasing 
the number of neutrons or protons, the effective spin-orbit interaction becomes weaker and 
this results in a reduction of the energy spacings for spin-orbit partners. The spin-orbit 
potential originates from the addition of two large fields: the field of the vector mesons 
(short range repulsion), and the scalar field of the sigma meson (intermediate attraction). 
In the first order approximation, and assuming spherical symmetry, the spin orbit term can 
be written as 

V s .o. = -^-V ls (r) , (28) 
r or 

where Vi s is the spin-orbit potential 

Tfi 

V ls = (V-S). (29) 

m ef.f 

V and S denote the repulsive vector and the attractive scalar potentials, respectively. m e // 
is the effective mass 

m eff = m- l -{V-S). (30) 

On the neutron-rich side the magnitude of the spin-orbit term V s , a . decreases as we add more 
neutrons, i.e. more units of isospin. This is reflected in the energy spacings between the 
neutron spin-orbit partner states 

A>Ei s = E nj ij = i-\/2 — E n jj = i+i/2 ■ (31) 

In Fig. H we plot the energy spacings between neutron spin-orbit partners in Sn isotopes, 
calculated in the RHB model with the DD-ME1 and NL3 effective interactions. The calcu- 
lated isotopic dependence is almost identical. Both interactions predict a reduction of the 
energy spacings between spin-orbit partners of ~ 50% in the interval 100 < A < 150. 
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In order to test the DD-ME1 effective interaction in the region of heavy nuclei, we have 
calculated the Pb isotopes with 196 < A < 214. In Fig. |1(] we display the deviations of 
the RHB theoretical masses of Pb isotopes from the empirical values fiOfl . The accuracy 
of the binding energies calculated with DD-ME1 is comparable to that obtained with the 
NL3 interaction. However, in contrast to the case of Sn isotopes, the DD-ME1 interaction 
systematically gives more binding as compared with NL3, especially for A < 208. 

Due to the intrinsic isospin dependence of the effective single-nucleon spin-orbit poten- 
tial, the relativistic mean-field models naturally reproduce the anomalous charge isotope 



shifts [§4J]. The well known example of the anomalous kink in the isotope shifts of Pb iso- 



topes is shown in Fig. |TT[ The results of RHB calculations with the DD-ME1 and NL3 
effective interactions, and with the Gogny D1S interaction in the pairing channel, are com- 



pared with experimental data from Ref. [45]. Both interactions reproduce the general trend 
of isotope shifts and the kink at 208 Pb. The effect is however, too strong with NL3. The 
experimental data are better described by the DD-ME1 interaction. 



Finally, in Fig. [12] we display the differences between radii of neutron and proton ground- 
state distributions of Pb isotopes, calculated with the DD-ME1 and NL3 effective interac- 
tions. Similar to the case of Sn isotopes, DD-ME1 systematically predicts much smaller 
values for r n — r p , in better agreement with available experimental data. The experimental 
values of r n — r p in 208 Pb are: 0.20 ± 0.04 fm deduced from proton scattering data |3"7|] , and 
0.19 ± 0.09 fm deduced from the alpha scattering excitation of the isovector giant dipole 
p8 |. In a recent analysis of intermediate energy nucleon elastic scattering data, 



resonance 



and correlated with analyses of electron scattering data, a value ~ 0.17 fm was suggested 
for the neutron skin thickness in 208 Pb [0. As it has been emphasized in a recent analysis 



of neutron radii in mean- field models [11], the value of r n — r p in 208 Pb is crucial for con- 
straining the isovector channel of effective interactions in the mean-field approach, both in 
non-relativistic and covariant models. 



V. PARITY- VIOLATING ELASTIC ELECTRON SCATTERING AND NEU- 
TRON DENSITY DISTRIBUTIONS 

Data on neutron radii and neutron density distributions provide not only basic nuclear 
structure information, but they also place additional constraints on effective interactions 
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used in nuclear models. Potentially, a very accurate experimental method for the deter- 
mination of neutron densities is the elastic scattering of longitudinally polarized electrons 
on nuclei. The parity-violating asymmetry parameter, defined as the difference between 
cross sections for the scattering of right- and left-handed longitudinally polarized electrons, 
produces direct information on the Fourier transform of the neutron density |41]]. A re- 



cent extensive analysis of possible parity-violating measurements of neutron densities, their 
theoretical interpretation, and applications can be found in Refs. ||| 



In Ref. [|5(| we have studied parity-violating elastic electron scattering on ground-state 
densities of neutron-rich nuclei that were calculated in the RHB model with the NL3 + 
Gogny D1S interaction. For the elastic scattering of 850 MeV electrons on these nuclei, the 
parity-violating asymmetry parameters were calculated using a relativistic optical model 
with inclusion of Coulomb distortion effects. The asymmetry parameters for chains of iso- 
topes were compared, and their relation to the Fourier transforms of neutron densities was 
studied. In this work we have shown that the new density- dependent effective interaction 
DD-ME1 predicts ground-state neutron density distributions that are in much better agree- 
ment with experimental data. Thus, in this section we include an analysis of parity- violating 
elastic electron scattering on 208 Pb and on those Sn isotopes for which there are data on 
r n — r p values. 

We consider elastic electron scattering on a spin-zero nucleus, i.e. on the potential 

V(r) = V(r)+ l5 A(r) , (32) 

where V(r) is the Coulomb potential, and A(r) results from the weak neutral current ampli- 
tude 

Mr) = §j- 2 Pw(r) • (33) 

The weak charge density is defined 

pw (r) = J d?r'G E {\v - v'\)[-p n {r') + (1 - Mn 2 e w )p p (r')} , (34) 

where p n and p p are point neutron and proton densities and the electric form factor of the 
proton is G_E;(r) « f^ e ~ Ar with A = 4.27 fin -1 , sin 2 Gv^ = 0.23 for the Weinberg angle. 
In the limit of vanishing electron mass, the electron spinor \l/ defines the helicity states 

tf ± = i(l±7 B )¥, (35) 
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which satisfy the Dirac equation 

[a-p+V ± (r)]V± = EV± , (36) 

with 

V±(r) = V(r)±A(r). (37) 

The parity- violating asymmetry Ai, or helicity asymmetry, is defined 

da + /dQ — da^/dQ 
1 = da + /dQ + da-/dQ ' (38) 

where +(— ) refers to the elastic scattering on the potential V±(r). This difference arises 
from the interference of one-photon and Z° exchange. 

Starting from the relativistic Hartree-Bogoliubov solutions for the self-consistent ground 
states, the charge and weak densities are calculated by folding the point proton and neutron 
densities. These densities define the Coulomb and weak potentials in the Dirac equation 
for the massless electron. The partial wave Dirac equation is solved with the inclusion of 
Coulomb distortion effects, and the cross sections for positive and negative helicity electron 
states are calculated. The parity-violating asymmetry parameters are plotted as functions 
of the scattering angle 6, or the momentum transfer q, and they are compared with the 
Fourier transforms of the neutron density distributions. 



4:71 f 

F(q) = — dr r 2 j (qr)p n (r) . (39) 
Q J 



In Fig. [13] we plot the parity-violating asymmetry parameters A[ for elastic electron 
scattering from 208 Pb at 850 MeV, as functions of the momentum transfer q = 2Esin6/2, 
and compare them with the squares of the Fourier transforms of the neutron densities 

Atx 
(1 

The solid and dotted curves correspond to RHB neutron ground state densities calculated 
with the NL3 and DD-ME1 effective Lagrangians, respectively. The asymmetries Ai are of 
order of < 10 -5 and increase with the momentum transfer q. We notice that, even though 
the values of r n — r p calculated with the two interactions differ by « 0.07 fm, the differences 
between the calculated asymmetry parameters Ai become more pronounced only for q > 1 
fm -1 . 

The parity-violating asymmetry parameters Ai for elastic scattering from even-A iso- 
topes 116 ~ 124 Sn at 850 MeV, as functions of the scattering angle 9, are shown in Fig. 14. 



They correspond to the ground-state densities calculated with the DD-ME1 interaction. In 
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Ref. [|5(| we have discussed the sensitivity of the asymmetry parameters to the formation 
of the neutron skin. For the Sn isotopes this effect is illustrated in Fig. |7], where the calcu- 
lated differences between neutron and proton radii of ground-state distributions are shown 
in comparison with experimental data. The corresponding asymmetry parameters < 10 -5 
in Fig. [TJ] display somewhat more pronounced differences between neighboring isotopes only 
for 9 > 15°. Finally, in Fig. [15] the calculated asymmetry parameters, as functions of the 
momentum transfer, are compared with the squares of the Fourier transforms of the neutron 
densities for 116_124 Sn. Obviously a resolution better than < for the asymmetry pa- 
rameters is necessary in order to obtain useful informations on neutron density distributions 
from parity-violating elastic electron scattering. This resolution might be already available 
at existing experimental facilities |fl9|| . For the Sn isotopes, in particular, the differences 
between neighboring isotopes are only seen for q > 1.5 fm _1 . We should also mention that 
the magnitude of the calculated asymmetry parameters depends, of course, on the electron 



energy. For electron energies below 500 MeV the asymmetry parameters are small |50 
while above 1 GeV the approximation of elastic scattering on continuous charge and weak 
densities is not valid any more, and the structure of individual nucleons becomes important. 

VI. SUMMARY AND CONCLUSIONS 

In the last couple of years the relativistic Hartree-Bogoliubov (RHB) model has been very 
successfully applied in the description of a variety of nuclear structure phenomena. With the 
standard non-linear meson-exchange relativistic mean-field effective interactions in the ph- 
channel, however, the predictive power of the RHB model is somewhat limited, especially 
for isovector properties of exotic nuclei far from /3-stability. We have tried to overcome 
these limitations by extending the RHB model to include density-dependent meson-nucleon 
couplings. The particular implementation of the model presented in this work is based 
on the density dependent relativistic hadron field (DDRH) theory |25| , |26|j . The effective 



Lagrangian in the p/i-channel is characterized by a density dependence of the a, uo and p 
meson-nucleon vertex functions. The single-nucleon Dirac equation includes the additional 
rearrangement self-energies that result from the variation of the vertex functionals with 
respect to the baryon field operators, and which are essential for the energy-momentum 
conservation and the thermodynamical consistency of the model. In this work we have 
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used the phenomenological density functional forms of the meson-nucleon coupling vertices 
introduced in the DDRH framework by Typel and Wolter |30|. The parameters of the new 
effective interaction DD-ME1 have been determined by a multiparameter fit constrained by 
properties of nuclear matter and by a set of experimental data on ground-state properties 
of spherical nuclei. Pairing correlations in the pp-channel of the RHB model are described 
by the pairing part of the finite range Gogny interaction. 

Properties of symmetric and asymmetric nuclear matter calculated with the new density- 
dependent effective interaction DD-ME1 have been compared to those obtained with the 



effective interaction of Typel and Wolter [p0| , and with two standard non-linear parameter 
sets sets NL3 [10| and NL1 It has been shown that the density dependent meson- 

nucleon couplings improve the behavior of the nuclear matter equation of state at higher 
densities and reproduce the empirical value of the asymmetry energy at saturation density. 
The properties of asymmetric nuclear matter are much better described by the two density 
dependent interactions and, in contrast to the non-linear NL3 and NL1 forces, these inter- 
actions reproduce the microscopic many-body neutron matter equation of state of Friedman 
and Pandharipande. 

The RHB model with the density-dependent interaction DD-ME1 in the p/i-channel, and 
with the finite range Gogny interaction D1S in the pp-channel, has been tested in the calcu- 
lation of ground-state properties of Sn and Pb isotopes. Results of fully self-consistent RHB 
calculations of binding energies, charge radii, differences between neutron and proton radii, 
spin-orbit splittings, performed with the interaction DD-ME1 and with the non-linear inter- 
action NL3, have been compared with available experimental data. While both interactions 
predict nuclear masses with the same level of accuracy (absolute deviations « 0.1 — 0.2%), 
the improved isovector properties of DD-ME1 result in a better description of charge radii, 
and especially the calculated values of r n — r p are in much better agreement with experi- 
mental data. The correct description of the data on differences of the radii of neutron and 
proton ground-state distributions, on neutron radii and neutron density distributions is very 
important for studies of new phenomena in exotic nuclei far from /3-stability (neutron skin, 
neutron halo, pygmy isovector dipole resonances), for astrophysical applications (properties 
of neutron stars, neutron capture rates), and for a theoretical interpretation of measure- 
ments of parity nonconservation effects in atomic systems (tests of the Standard model of 
electroweak interactions) j5T|. In principle, very accurate data on neutron density distri- 
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butions could be obtained from the elastic scattering of longitudinally polarized electrons 
on nuclei. Using the ground-state densities calculated with the DD-ME1 interaction in the 
RHB model, we have performed an analysis of parity-violating elastic electron scattering 
on 208 Pb and on those Sn isotopes for which there are experimental data on r n — r p . For 
the elastic scattering of 850 MeV electrons on these nuclei, the parity-violating asymmetry 
parameters have been calculated using a relativistic optical model with inclusion of Coulomb 
distortion effects, and related to the Fourier transforms of the neutron density distributions. 

The RHB model with density-dependent meson-nucleon couplings represents a significant 
improvement in the relativistic mean-field description of the nuclear many-body problem 
and, in particular, of exotic nuclei far from /3-stability. The improved isovector properties 
of the effective interaction in the p/i-channel on one hand, and the unified description of 
mean-field and pairing correlations in the Hartree-Bogoliubov framework on the other, offer 
a unique possibility for accurate studies of nuclei with extreme ground-state isospin values 
and with Fermi levels close to the particle continuum. Particularly interesting will be studies 
of deformed nuclei for which unusual shape coexistence phenomena, and even very different 
proton and neutron ground-state deformations, are expected far from stability and close to 
the drip-lines. Isovector ground-state deformations could also give rise to exotic modes of 
low-energy isovector collective excitations. We have already started with density-dependent 
RHB calculations of deformed nuclei, and work is also in progress on the description of 
collective excitations in the framework of relativistic RPA/QRPA with density-dependent 
interactions. 
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TABLE I: The effective interaction DD-ME1. The masses and meson- nucleon couplings are shown 
in comparison with the parameters of the density-dependent mean- field model of Ref. (3(| (TW-99). 





UJJ-MiIjI 


1 W-99 




549.5255 


550.0000 


m u 


783.0000 


783.0000 


m p 


763.0000 


763.0000 


ga(Psat) 


10.4434 


10.7285 


9uj(Psat) 


12.8939 


13.2902 


9p{Psat) 


3.8053 


3.6610 


a a 


1.3854 


1.3655 


K 


0.9781 


0.2261 




1.5342 


0.4097 


da 


0.4661 


0.9020 


aui 


1.3879 


1.4025 




0.8525 


0.1726 




1.3566 


0.3443 




0.4957 


0.9840 


a p 


0.5008 


0.515 



TABLE II: Nuclear matter properties calculated with the density-dependent effective interactions 
DD-ME1 and TW-99 and the non-linear parameter sets NL3 @ and NL1 ||. 





DD-ME1 TW-99 


NL3 


NL1 


Psat (fm 3 ) 


0.152 


0.153 


0.149 


0.153 


E/A (MeV) 


-16.20 


-16.25 


-16.25 


-16.42 


K (MeV) 


244.5 


240.0 


271.8 


211.3 


m* 


0.578 


0.556 


0.60 


0.57 


a 4 (MeV) 


33.1 


32.5 


37.9 


43.7 


po (MeV/fm 3 ) 


3.26 


3.22 


5.92 


7.0 


AK (MeV) 


-128.5 


-126.5 


52.1 


67.3 
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TABLE III: Binding energies, charge radii and differences between neutron and proton radii used 
to adjust the parameters of the DD-ME1 interaction. The calculated values are compared with 
experimental data (in parentheses). 







^ivre v ) 


r c h (fm) 


r n — r p (fm) 


lbQ 


-7.974 


(-7.976) 


2.730 (2.730) 


-0.03 


40 Ca 


-8.576 


(-8.551) 


3.464 (3.485) 


-0.05 


48 Ca 


-8.631 


(-8.667) 


3.482 (3.484) 


0.19 


9U Zr 


-8.704 


(-8.710) 


4.294 (4.272) 


0.06 


1 1 O 

112 Sn 


-8.501 


(-8.514) 


4.586 (4.596) 


0.11 


116 Sn 


-8.516 


(-8.523) 


4.616 (4.626) 


0.15 (0.12) 


124 Sn 


-8.462 


(-8.467) 


4.671 (4.674) 


0.25 (0.19) 


132 Sn 


-8.352 


(-8.355) 


4.720 


0.27 


204p b 


-7.885 


(-7.880) 


5.500 (5.486) 


0.18 


208p b 


-7.884 


(-7.868) 


5.518 (5.505) 


0.20 (0.20) 


214p b 


-7.764 


(-7.772) 


5.568 (5.562) 


0.27 


210p o 


-7.857 


(-7.834) 


5.553 


0.18 
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TABLE IV: Energy separation (in MeV) between spin-orbit partner states in doubly closed-shell nu- 
clei, calculated with the DD-MEl and NL3 interactions, and compared with experimental data [39 J. 







DD-MEl 


ATT O 

JNL3 


Exp. 


16q 


vlp 


6.316 


6.482 


6.18 




■nip 


6.249 


6.404 


6.32 


40 Ca 


vld 


6.567 


6.716 


6.00 




nld 


6.507 


6.630 


6.00 


48 Ca 


vlf 


7.689 


7.542 


8.38 




v2d 


1.723 


0.888 


2.02 


132 Sn 


v2d 


1.883 


1.573 


1.65 




nig 


6.244 


6.230 


6.08 




n2d 


1.822 


1.584 


1.75 


208p b 


1/2/ 


2.197 


1.860 


1.77 




vli 


6.839 


6.813 


5.84 




v3p 


0.878 


0.802 


0.90 




n2d 


1.647 


1.525 


1.33 




TTlh 


5.837 


5.809 


5.56 
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FIG. 1: Density dependence of the couplings of the a-, cj-, and p-meson. The result of the 
present analysis (DD-ME1) is shown in comparison with the parameters of the effective interaction 
TW-99 @. 
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FIG. 2: Binding energy per nucleon for symmetric nuclear matter as a function of the baryon 
density, calculated with the density-dependent effective interactions DD-ME1 and TW-99 [ppf] , and 
the non-linear parameter sets NL3 pi and NL1 p3f. 
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FIG. 3: S2 (p) coefficient ( p7| ) of the quadratic term of the energy per particle of asymmetric 
nuclear matter, calculated with the effective interactions DD-ME1, TW-99 [30|, NL3 [1C], and 
NL1 @. 
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FIG. 4: Energy per particle of neutron matter as a function of the neutron density. The results 



for the four relativistic mean-field interactions DD-ME1, TW-99 J30|, NL3 g(J, and NL1 J33| are 
shown in comparison with the neutron matter equation of state of Friedman and Pandharipande 
(FP-81) C 
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FIG. 5: The deviations (in percent) of the theoretical masses of Sn isotopes, calculated in the RHB 



model with the DD-ME1 and NL3 interactions, from the empirical values [40|. The pairing part 
of the Gogny interaction D1S has been used in the pp channel of the RHB model. 
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FIG. 6: The deviations (in percent) of the theoretical charge radii of Sn isotopes, calculated in 



the RHB model with the DD-ME1 and NL3 interactions, from the experimental values [41]. 
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FIG. 7: Differences between neutron and proton radii of ground-state distributions of Sn isotopes, 
calculated with the DD-ME1 and NL3 effective interactions. 
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FIG. 9: Energy spacings between neutron spin-orbit partner states in Sn isotopes, calculated in 
the RHB model with the DD-ME1 and NL3 effective interactions. 
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FIG. 10: The deviations (in percent) of the theoretical masses of Pb isotopes, calculated in the 



RHB model with the DD-ME1 and NL3 interactions, from the empirical values [4C]. 
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FIG. 11: Charge isotope shifts in even-A Pb isotopes. The results of RHB calculations with 
the DD-ME1 and NL3 effective interactions, and with the Gogny D1S interaction in the pairing 



channel, are compared with experimental data from Ref. [45]. 
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FIG. 12: Differences between neutron and proton radii of ground-state distributions of Pb isotopes, 
calculated with the DD-ME1 and NL3 effective interactions. 
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FIG. 13: Parity-violating asymmetry parameters A\ (upper panel) and squares of normalized 
Fourier transforms of neutron densities (lower panel), as functions of the momentum transfer q, for 
elastic scattering from 208 Pb at 850 MeV. 



40 




e[deg] 

FIG. 14: Parity-violating asymmetry parameters A[ for elastic scattering from 116 ~ 124 Sn at 850 
MeV, as functions of the scattering angle 9. 
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FIG. 15: Parity-violating asymmetry parameters A\ (upper panel) and squares of normalized 
Fourier transforms of neutron densities (lower panel), as functions of the momentum transfer q, for 
elastic scattering from 116 ~ 124 Sn at 850 MeV. 
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